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(57) Abstract: The present invention relates to an installation for preventive identification of earthquakes, composed of multiple 
peripheral stations capable of detecting the very low frequency electromagnetic waves originated because of the mechanical defor- 
mation of underground layers of rocks and sending results to a central unit that processes received data in order to obtain the epicentre 
coordinates and depth of a future earthquake. 
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Description 

Installation for preventive identification of earthquake epicentre and 
depth. 

The present patent application relates to an installation for preventive 
identification of earthquake epicentre and depth based on the detection and 
processing of the very low frequency (VLF) electromagnetic wave emissions 
which forewarn earthquakes, 
5 The invention, of which no antecedents are known, is based on the 

fact that in earthquake zones the earth's crust is periodically subject to strong 
tensions and compressions that originate earthquakes when they reach the 
breaking point. 

It is known that the geological structure of earthquake zones of the 
10 earth's crust is mainly composed of rocks with high contents of crystals. 
When subject to strong mechanical deformation, rocks produce electricity 
(piezoelectricity) with consequent generation of VLF electromagnetic waves, 
whose intensity and frequency are related to the mechanical tensions that 
originated them. 

15 The frequencies of the electromagnetic waves range from fractions of 

Hertz to some tens of kilohertz; since the soil tends to attenuate higher 
frequencies, the deeper the phenomenon that originated the waves, the 
lower the received spectrum will be. 

Being known that the intensity of electromagnetic emissions is 
20 proportional to the square of the forces applied by volume unit of the 
piezoelectric rock, an accurate statistical study of the territory will allow to 
predict the intensity of the future earthquake with good approximation. 

It must be noted that the electromagnetic waves generated in the 
mechanically anomalous area attenuate with the square of the distance 
25 between the origin and reception point. 

This means that the farther we get from the source, the weaker the 
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intensity of forewarning emissions will get and the same will be hidden by 
man-generated interference, static discharges and foreign radio waves. 

The main purpose of the present invention is to predict the location of 
the earthquake epicentre and depth based on the analysis of forewarning 
■5 VLF electromagnetic waves. 

The earthquake epicentre is determined by using a network of 
peripheral receiving stations (minimum three) located on the territory at a 
distance of few hundreds of kilometres, possibly far away from interference 
sources, such as high voltage electric lines, towns and industrial centres. 
10 Each peripheral station is composed of a series of suitable 

components (as illustrated with details in the description below) housed in a 
metal container. 

Peripheral stations are programmed to automatically send results to a 
central unit over radio or telephone. The centra) unit makes a special 
15 innovative comparison (as described below) of the signals received from the 
peripheral stations in order to locate the epicentre of the future earthquake. 

For major clarity the description of the invention continues with 
reference to the enclosed drawings, which are intended for purposes of 
illustration and not in a limiting sense, whereby: 
20 - Fig. 1 is a diagrammatic view of a peripheral station used to detect VLF 
electromagnetic waves; 

- Fig. 2 is a block diagram of the operating components of the peripheral 
station; 

- Fig. 3 is a diagram that illustrates the method used by the central 
25 monitoring unit to process the information received from, the peripheral 

stations in order to locate the earthquake epicentre; 

- Fig. 4 is a three-dimensional view of a portion of "the earth's crust with four 
peripheral stations according to the installation of the invention. 

With reference to Figures 1 and 2, each peripheral station makes use 
30 of two receiving antennas: a main antenna (1 ) with balanced screened loop 
and a preferably underground service wire antenna (2). 
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The signals emitted from the two antennas (1, 2) are conveyed 
towards suitable devices housed in the container (C), after being adapted in 
impedance, amplified and filtered with suitable filters (3) in order to eliminate 
interference deriving from external interference sources, 
5 Once inside the container, the signals detected by the antennas (1, 2) 

are converted from analogue into digital by means of suitable A/D devices 

(4) . 

The signals, now in digital format, are sent to a processing computer 

(5) after being marked with a clock signal (6) of one PPS (pulse per second) 
10 coming from a GPS satellite receiver (7) that supplies sample-time. 

The computer (5) starts digital signal processing (DSP) in order to 
suppress the interference of the VLF electromagnetic waves through the 
comparison of the signals received by the two antennas (1 , 2). 

The signal obtained at the end of the digital signal processing (DSP) is 
15 memorised in a suitable memory (8) and conditioned for the subsequent 
transmission to a central monitoring and processing unit - not shown in the 
enclosed figures - over a telephone line connected to a modem (9a) or radio 
modem (9b) associated with a transmitting antenna (A) mounted on the roof 
of the container (C), 

20 As soon as the attention level is exceeded by a peripheral station, the 

central monitoring and processing unit acquires the data from the stations via 
radio or telephone in real time and starts processing. 

At this regard, the central unit is equipped with a first computer used to 
acquire and memorise the signals coming from the stations located on the 
25 monitored territory, and a second computer used for DSP (Digital Signal 
Processing) operations in order to process the signals received from the 
peripheral stations. 

The operating logic of the second computer is diagrammatically 
illustrated in Figure 3, which shows the processing of signals coming from 
30 four peripherat stations (Rx-A, Rx-B, Rx-C, and Rx-D). 

More precisely, the operation consists in processing the signals 
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coming from the peripheral stations, after suitable equalisation, to find the 
correlation peak (the correlation consists in the temporal flow of the signals 
received from the peripheral' stations in numerical mode until they coincide 
within a predefined time period). 
'5 The correlation peak is used for the temporal alignment of the signals 

received from the peripheral stations in order to measure the relative delays 
ft(a-b), t(a-c), t(a-d)] between the time marks of the signals received from the 
stations located in the territory. 

In other words, the VLF signals received from the stations are marked 
10 with a time signal derived from a GPS, and time differences of arrival are 
calculated in the data collection and processing centre with the correlation 
and GPS marking technique, as shown in Figure 3. 

Delays are measured on several signal samples received from each 
station. The results are then averaged and data is acquired when minimum 
1 5 variance is achieved. 

Now, by knowing the station coordinates and the signal time difference 
of arrival, the epicentre coordinates and the earthquake depth can be 
obtained by using the TDOA (Time Difference Of Arrival) technique. 

This technique allows for the 2-D identification (x-y) of the earthquake 
20 epicentre with three stations, in addition to higher accuracy, the use of four or 
more stations allows for the possibility of 3-D operation, which includes the z 
axis representing the epicentre depth. 

With reference to a Cartesian system (x, y, z) as shown in Fig. 4, the 

■ 

unknown coordinates of the emission signal are obtained by solving the 
25 system: 



30 



V(*-xJ + (y-y b f + {z~z b f -4{x-x a f + (y~yj +(z-z a f =c^ 

+Cv-^) 2 +(*-0 2 -V(*-*J 2 +(y-yJ +(*-0 2 =ct ae 
+(y-y*) 2 Az-zrf +b-yJ =c-^ 
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where t(ab), t(ac) t t(ad) are the delays of the signals received from stations B, 
C and D, respectively, with respect to the signal received from station A and 
constant C is the propagation speed of electromagnetic waves underground. 
Therefore, the Cartesian system can be used to identify the epicentre 

^5 coordinates (x, y), as well as its depth (2). 

In case of three peripheral stations, the processing procedure should 
consider the delays of the signals from two stations with respect to the 
signals from the third station, and the calculation system should be 
composed of only two equations with two unknowns, that is to say the 

10 epicentre co-ordinates (x, y). 
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Claims 

1) Installation for preventive identification of earthquakes, characterised in 
that it is composed of multiple peripheral stations capable of detecting the 
very low frequency (VLF) electromagnetic waves originated because of the 
mechanical deformation of the underground layers of rocks and sending 

-5 results to a central processing unit; it being provided that each peripheral 
station is composed of two receiving antennas: a main antenna with balanced 
screened loop (1) and a preferably underground service wire antenna (2), 
whose signals, after being adapted in impedance, amplified and filtered, are 
converted from analogue to digital by means of suitable A/D devices (4), and 

10 then sent to a processing computer (5) after being marked with a clock signal 
(6) of one PPS (pulse per second) coming from a GPS satellite receiver (7) 
that gives the sample-time; it being provided that the computer (5) starts 
digital signal processing (DSP) for the suppression of interference in VLF 
electromagnetic waves thanks to the comparison of the signals received from 

15 the two antennas (1, 2), so that the processed signal can be memorised in a 
memory (8) and sent to the central processing unit via telephone or radio as 
soon as the attention level of the peripheral station is exceeded; it being 
provided that the centra! processing unit is equipped with a first computer 
used to acquire and memorise the signals (Rx-A, Rx-B, Rx-C, Rx-D Rx-n) 

20 coming from the stations located in the monitored territory, and a second 
computer used to process the signals, after equalisation, to find the 
correlation peak according to the method illustrated in the description above, 
which allows for the temporal alignment of the signals in order to measure 
the relative delays [t(a-b), t(a-c), t(a-d)„..t(a-n)] between the time marks of 

25 the signals received from the peripheral stations;' it being provided that the 
central processing unit is capable of measuring the delays on multiple signal 
samples received from each peripheral station and averaging the results so 
that data is acquired and used when minimum variance is achieved, knowing 
the station coordinates and the signal time differences of arrival, in order to 
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obtain the epicentre co-ordinates and the depth of the future earthquake, 
based on the processing technique known as TDOA (Time Difference Of 
Arrival) and on the calculation modes illustrated in the description above. 
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